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We present a general procedure for applying the scale-setting prescription of Brodsky, Lepage and Mackenzie to 
higher orders in the strong coupling constant a s . In particular, we show how to apply this prescription when the 
leading coefficient or coefficients in a series in a s are anomalously small. We give a general method for computing 
an optimum scale numerically, within dimensional regularization, and in cases when the coefficients of a series are 
known. We find significant corrections to the scales for R e + e - , F(B — > X u eV), F(t — » bW), and the ratios of the 
quark pole to MS and lattice bare masses. 



1. Introduction 

QCD processes computed to a finite order in 
perturbation theory depend on the scale chosen 
for the running coupling constant a s (q). While 
these variations diminish as higher orders are in- 
cluded, for low-order calculations they can be sig- 
nificant. Finding an optimum, physically moti- 
vated method for choosing this scale is important 
not only to produce accurate results, but also to 
reasonably estimate convergence based on the size 
of the series coefficients. Such a method allows a 
meaningful prediction or comparison with data 
even at leading order. 

In this paper, we apply the prescription de- 
fined by Brodsky, Lepage and Mackenzie (BLM) 
in Refs. [jjJH beyond leading order. In so doing, 
we remedy an anomaly observed in a variety of 
applications as discussed, for example, in Ref. ||. 
We show that this requires a simple extension of 
the calculation needed to set the scale at lowest 
order. We find that it leads to significant correc- 
tions in the scales of several important processes, 
and allows us to extract masses from lattice sim- 
ulations which were previously inaccessible. (See 
Refs. for other higher-order extensions.) 
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2. The Prescription 
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Figure 1. The BLM prescription fixing the opti- 
mum scale q* to leading order in a s (<?*). Vacuum 
bubbles include both fermions and gluons. 



The original statement of the prescription, 

«.(«*) / d 4 qf(q) « J d 4 qa s (q)f(q) W 

illustrated in Fig. |l|, is to choose the scale q* 
associated with a particular gluon such that it 
matches, as well as possible, the same diagram 
with a fully dressed gluon [Q. Here f(q) is the 
integrand within the diagram which includes the 
gluon. In this way, a s (q*) absorbs much of the 
effect of these higher-order corrections, and rep- 
resents the true strength of this gluon's coupling 
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to the rest of the diagram. When f(q) is sensitive 
to large q, expanding a s on both sides at a com- 
mon scale (i gives a form suitable for numerical 
calculation, 

M,-^)= (/fa) ;°y 2)> (2> 

s «log(, 2 /f J ))) , 

with q* given by the average log of the momentum 
running through that gluon [g| . 

3. The Problem 

When (f{q)) vanishes, Eq. (||) is meaningless. 
Even when nonzero, if (f{q)) is anomalously small 
compared to (/(g) log(q 2 /^i 2 )), the result can be 
misleading. This occurs in a variety of processes, 
and for some set of parameters in most processes. 
For example, it occurs for the perturbative re- 
lation between the bare lattice quark mass in 
NRQCD and its pole mass at the value rele- 
vant for charm. A gluonic contribution for which 
(f{q)) — ► is dominated by its second-order term; 
Eq. (g) fails because it attempts to match a first- 
order statement of the general prescription to a 
process which is properly second-order. 

4. The Solution 

The solution is to apply the prescription be- 
yond leading order, as illustrated in Fig. [| This 
gives a reasonable q* when the leading diagram is 
anomalously small, and it approaches a sensible 
limit when it vanishes altogether. In rare cases 
when both the first and second order contribu- 
tions are small, it is simple to carry this one loop 
higher. 

The result of matching at a common scale /j. 

is a 
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log(g* 2 /M 2 ) = ((log^ 2 /^ 2 ))) ± 
with 

CT 2 EE «k,gV/ M 2 ))) - ((log^ 2 /^ 2 ))) 2 . (4) 

When f(q) does not change sign, it may be inter- 
preted as a probability distribution, so a 2 > 0, 
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Figure 2. The BLM prescription applied to sec- 
ond order in a s (q*). 



and Eq. (g) gives the correct scale. However, 
when a 2 < 0, there are evidently significant sign 
changes, and (/) will be anomalously small. In 
this case, Eq. (|^) gives the correct choice of scale. 
It causes the leading-order discrepancy between 
the left and right sides in Fig. ^ to vanish, and 
minimizes it at next order Applying Eq. (^) 
requires only the additional computation of the 
average log squared within the same integrand as 
in Eq. (|). 

5. A Simple Example 

Consider as a simple model for the integrand 

f(q) = (l + c)S(q-q a )-S(q-q b ), (5) 

sensitive to two scales, q a and q^. It suffers from 
partial cancellations when c > —1, and a 2 be- 
comes negative; when c = 0, (/) vanishes iden- 
tically and Eq. (j^) blows up. In Fig. [| we show 

and 
be- 
haves as expected, with q* approaching q a when 
|c| 3> 1, qb when c = —1, remaining between in 
the interim. Even in the region where c is large 
and positive and Eq. (||) gives a well-behaved q* , 
Eq. (0) is clearly preferable. 
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the result of applying Eq. (H) in this region, ; 
Eq. (||) when c < — 1 and a 2 > 0. It clearly 



6. Determining q* in MS 

The g4> 3 self-energy diagram illustrated in 
Fig. H presents a more realistic example. Evaluat- 
ing the diagram with an additional denominator 
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(2tt)™ k 2 +m 2 (p-k) 2 +m 2 (fc 2 //z 2 )' 3 



(6) 
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Figure 3. The BLM scale q* as a function of c, 
with q a = 2.0 and qb = 1.8. The first order solu- 
tion determines q* for c < — 1, second order for 
c > — 1. Dark dotted lines show the first-order so- 
lution in regions in which it does not apply; light 
dotted lines display inapplicable second-order so- 
lutions. 



and expanding the result in (3 produces the se- 
quence of average logs needed to set the scale to 
any desired order. This is not much more diffi- 
cult within dimensional regularization than com- 
puting the diagram itself. In Fig. || we display 
the result for q* as a function of the propagator 
momentum p Q. The second-order solution is 
appropriate in the intermediate region. It gives a 
physically reasonable result which connects con- 
tinuously with the first-order solution appropriate 
for large and small p/m, whereas the first-order 
solution diverges in this region. 




P P 



Figure 4. One loop self-energy diagram in the 
gcj) 3 model. 



Figure 5. The MS BLM scale q*/p as a function 
of momentum p/m for the 4> 3 self-energy diagram. 



7. Determining q* from an Existing Series 

When the terms in a series are already known, 
the average logs needed for scale setting can be 
extracted by using the dependence on the num- 
ber of flavors n/, or equivalently, /3q = (11 — 
(2/3)n/)/47r. Contributions from a particular 
dressed gluon within a diagram take the form 

c a s (/x) + (oi - ci/3 )a%(ji) 

+ (oa + • • • + c^l)al{n) + ■ ■ ■ . (7) 

The highest power n of /3q at each order in a s is 
due to n one- loop vacuum polarization diagrams, 
and so at large q we may make the identification 

co = (/) (8) 
cx/cq * «log(g 2 /M 2 ))) 
c 2 / C0 « ((log 2 (<z 2 /^ 2 ))) , 

allowing us to apply second-order scale setting 
for a series known to af. Below we use this to 
present several processes for which Eq. (H|) gives 
the appropriate scale. 
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8. Higher orders 

The series coefficients from a single dressed 
gluon depend on the scale q* according to 0] 

+ a 3 s (q*)[P 2 A 2 +(3 1 A 1 ] 
+ ai(q*)[f A 3 + ^lAa + ^Ax] 

+ a 5 a (q*MA4+™p 2 (3 1 A 3 

+ 3/3 /3 2 A 2 + ^I3\A 2 + /%A X ] + •••}, (9) 

with 

A„ ee (( [log^* 2 /^) - log(g 2 /M 2 )] "» • (10) 

In this form, the objective in choosing an opti- 
mum scale is transparent, regardless of the num- 
ber of loops kept in a s . For a gluon sensitive 
to a narrow, large momentum region, choosing 
the correct q* will approximately minimize the 
moments in Eq. (|i"o|). Contributions from the 
higher order diagrams which dress that gluon will 
be small, having been largely absorbed into the 
leading term. However, when the region of sensi- 
tivity in momentum is large, these moments will 
be minimized less effectively by a single scale q*, 
and higher order contributions will be larger. 

The general prescription to any order is then 
to skip any leading terms which are anomalously 
small, which will rarely be other than the first; 
choose the q* which eliminates the coefficient af- 
ter the first nonanomalous term and minimizes 
the next; and use higher moments when available 
to check the consistency of the scale choice. 

In Fig. H we display the dependence of these 
moments on log(g* 2 ) for the semileptonic de- 
cay width T(B — > X u ev e ) expressed in terms 
of the MS b mass, using results from Ref. J?]]. 
In this case, Eq. (^) gives the appropriate scale, 
and examination of the higher moments con- 
firms this. Choosing the second-order prescrip- 
tion eliminates the second moment, minimizes the 
third moment, and log(g* 2 ) is near either the zero 
or minimum for all the higher moments. As a re- 
sult, higher order coefficients are small, and the 



leading term effectively represents the strength of 
this gluon's coupling. 
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Figure 6. The moments lAnl^™ as functions of 
log(g*7^ 2 ) for n = 1 to 8 for T(B -> X u eV e ). 
The vertical line indicates the preferred scale from 
Eq. (|). 



9. Applications 

When at least the second logarithmic moment 
is available, we may apply scale setting beyond 
leading order. Ref. [Q presents a collection of such 
series. In Table 111 we list results for four of these 
for which Eq. (0) gives the appropriate scale: 
R e + e -, the ratio of the pole to MS mass M/M, 
and the B and t decay widths T(B — ► X u ev) and 
r(i — > bW) expressed in terms of MS masses. 
These new scales represent significant corrections 
to the scales set by Eq. (||). While the new scale 
for M/M is increased, we note that the range of 
important momenta Aq is still relatively large, in- 
dicating sensitivity to low-momentum scales even 
when M is large, and threatening growing higher- 
order coefficients. This is not the case for B and 
t decays provided their series are expressed in 
terms of the short-distance MS masses, as indi- 
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Table 1 

Applications of second-order scale setting to several processes. q% gives the scale set by Eq. (||); q2 gives 
the preferred scale by Eq. (H). Aq, measures the range of momentum running through the gluon|4|. 



ci/co 



91 



C2/C0 



92 



A? 



-.691772 0.7076^ 

M/M: 

-4.6862 0.09603M 

£(g -> X u eu)/M 5 b : 

-4.3163 0.11554M 6 
r(i -► ftWO/ftiff : 



-5.7076 



0.05763M, 



-.186421 
17.623 
8.0992 
6.0996 



-.66497 
-4.3374 

-10.531 
-26.477 



1.06Vs 
0.27205M 

0.58534M 6 

0.75502M t 



0.38M 
0.35Mb 
0.34M, 



as 



cated by their ranges Aq and their moments 
in Fig. |. 

In Fig. |?j, we show the scale associated with 
the series connecting the quark pole mass to the 
NRQCD lattice bare mass. At both small and 
large values of the bare mass, Eq. (||) is appropri- 
ate and necessary to give an accurate measure of 
the optimum scale. 




Figure 7. The BLM scale q* for the pole mass 
as a function of bare lattice mass 0M0. The first 
order solution determines q* between aMo = 2.00 
and 3.50, second order elsewhere. Dotted lines 
show the first-order solution in regions in which 
it does not apply. The point at aMo = 7.00 is off 
the plot. 



10. Conclusions 

We have presented a simple general procedure 
for applying BLM scale-setting beyond leading 
order. Our main result is the second-order pre- 
scription Eq. (||), which is appropriate when the 
leading contribution is anomalously small. We 
gave a general method for computing an optimum 
scale numerically, within dimensional regulariza- 
tion, and in cases when the coefficients of a series 
are known, and applied it to several processes. 
Finally, we discussed its application at higher or- 
ders. 
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